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Abstract 

Heat resources of small temperature difference are easily accessible, free and 
unlimited on earth. Thermoelectric effects provide the technology for con- 
verting these heat resources directly into electricity. We present designs 
of electricity generators based on thermoelectric effects and using heat re- 
sources of small temperature difference, e.g., ocean water at different depths 
and geothermal sources, and conclude that large-scale power plants based 
on thermoelectric effects are feasible and economically competitive. The key 
observation is that the power factor of thermoelectric materials, unlike the 
figure of merit, can be improved by orders of magnitude upon laminating 
good conductors and good thermoelectric materials. The predicted large- 
scale power plants based on thermoelectric effects, if validated, will have a 
global economic and social impact for its scalability, and the renewability, 
free and unlimited supply of heat resources of small temperature difference 
on earth. 

Keywords: Thermoelectric Power Plants, Power Factor, Ocean and 
Geothermal Energy 



1. Introduction 

By the Second Law of thermodynamics, heat, as a form of energy, can 
be converted into usable energy only if there exists a temperature different 
between two heat reservoirs. The ideal efficiency of conversion is given by 
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the Carnot efficiency: 



^carnot = ™ j \*-) 
i-h 

where and T c are the high and low temperature, respectively. From this 
viewpoint, the largest accessible energy reservoir on earth is undoubtedly the 
oceans in regard of the thermocline in oceans: the surface water temperature 
within 100 meters of sea level is around 20K higher than that below 600 
meters in tropical regions since surface water is heated by, stores and con- 
centrates solar energy If this energy resource can be converted to electricity 
at a competitive cost, even with a fraction of the ideal Carnot efficiency, the 
energy problem in the world is solved once and forever. 
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Figure 1: A thermoelectric power plant using temperature difference of ocean water 
at different depths: cold water is pumped and passes through the heat exchanger 
/ generator where it is heated by surface water and generates electricity. Energy 
of ocean waves can be easily harvested to drive the flows. Surface water has 
higher temperature since it absorbs, stores and concentrates sunlight energy. The 
temperature profile is sketched according to pQ in tropical regions. 

Conventional technologies of large-scale power plants burning chemical 
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or nuclear fuels have to first convert heat into kinetic energy, which is then 
used to propel electromagnetic power generators. These power plants cannot 
make use of heat resources of small temperature difference. Thermoelectric 
(TE) materials, on the other hand, can directly convert heat into electricity 
no matter how small the temperature difference is [2] . Small-scale TE genera- 
tors have been commercialized in applications of microelectronics for reliable 
and independent power [3l H] , automobiles for improving overall fuel efficiency 
[5], and electricity supply in remote areas [6]. In conventional designs, the 
efficiency of conversion is a predominant design criterion to control cost since 
fuels account for the major cost of power generation [TJ El El QUI El [12] • The 
efficiency of a TE generator is dictated by a material property, namely, the di- 
mensionless figure of merit which, to our best knowledge, is hard to improve, 
and the state-of-the-art TE materials have a dimensionless figure of merit of 
around 2.4 at the room temperature [131 E]> corresponding to an efficiency 
of about 20% of the ideal Carnot efficiency. A larger dimensionless figure 
of merit is necessary for large-scale TE generators to be competitive against 
conventional conversion methods if fuels are consumed, which motivates re- 
cent renewed interest in thermoelectric materials focusing on improving the 
conversion efficiency. 

Here we present designs of TE generators at the scale of megawatt. The 
heat resources are taken as free and unlimited, e.g., ocean water at different 
depths (Fig. [TJ |1J), geothermal sources (Fig. |2j [2]), or waste heat from 
conventional power plants [16J. The main cost of the generator comes from 
the materials, the fabrication and assembly processes of the generator, and 
possibly the energy needed to drive the flows for heat exchange and electric- 
ity generation. This cost basis, as will be shown below, is realistic and shifts 
the importance of conversion efficiency (figure of merit) to conversion capac- 
ity (power factor), and favors simple structures for ease of mass production. 
Power factors of thermoelectric materials, unlike figure of merit, can be im- 
proved by orders of magnitude by engineering thermoelectric structures, e.g., 
laminating two or more TE materials. 

Below we focus our analysis on large-scale ocean-based TE power plants 
using water at different depths as "fuel" for its renewability and scalability. 
As illustrated in Fig. [TJ cold water is pumped and driven through tubes 
in a heat exchanger where it is heated by surface water. The tubes are 
made of TE materials, and therefore, as heat exchanges across the TE tube 
walls, electricity is generated. In spite of the well-known shortcoming of 
low conversion efficiency, large-scale ocean-based TE power plants enjoy the 
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following unique and unprecedented advantages: 



(i) Heat resources of small temperature difference, i.e., water at deep sea 
(below 600m), are easily accessible, free and unlimited. The power 
cost to pump water and drive the flows can be made free by harvesting 
energy from ocean waves. A valve can be designed to open in a valley 
(below sea level) and close on a peak (above sea level) of ocean waves, 
and henceforth an average water height above sea level in the cold water 
tank can be maintained without any power cost. 

(ii) It costs no land and can be installed in three dimensions, and hence 
has superior scalability. 

(iii) It is green, reliable since there is no moving solid parts, and has low 
maintenance cost. 

(iv) The power output is independent of the hour of a day and season at a 
tropical site. 

(v) It counteracts global warming. 

In addition we notice that ocean waves and convection naturally mix 
surface water and hence no power is needed to drive the flows exterior to 
the heat exchanger, and that there is no stringent structural requirement 
on the long underwater pipe if the material density is close to that of water. 
Structural materials of this kind can be realized by reinforced foam or plastics 
at a low cost which also serve as good thermal insulators. 

2. Proposed Designs of the generator 

Below we show large-scale ocean-based TE power plants at the order of 
megawatt are feasible and economically competitive even if power is needed 
to pump water and drive the flows on a waveless day. To see this and for 
simplicity, we assume the flows in the underwater pipe and TE tubes are 
laminar, and classic fluid dynamics asserts that the power needed to pump 
water to a height h and the flow rate through TE tubes driven by a pressure 
pgh is given by [17j: 
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Figure 2: A thermoelectric generator using geothermal resources. Hot water is 
pumped and passes through the heat exchanger where it is cooled by air and 
generates electricity. 



where where /i = 1.0 x Pa-s (p = 10 3 kg/m 3 ) is the viscosity (density) 
of water, N pp (iV tb ) is the number of underwater long pipes (near surface 
short TE tubes), L pp and A pp (L t b , Ab) are the length and cross-sectional 
area of the pipe (TE tube), Q pp (Qtb) is the volumetric flow rate through the 
pipe (TE tube), and 7 PP (7tb) is a dimensionless geometric factor depending 
on the cross-sectional shape of the pipe (TE tube). If both are circular, 



7 PP = 7tb 



8tt. Further, N pp Q 



pp^pp 



N th Q 



tb 



Q is the total volumetric 



flow rate through the generator, the second term in Pi oss can be identified 
as pghQ, and Pi oss scales as l/N pp and 1/iVtb for fixed Q. The formula (j2]) 
underestimates the actual power loss for neglecting turbulence, friction with 
the walls, etc, but is accurate at the order of magnitude and sufficient for 
our subsequent analysis. 

Electrically, the TE tube shall be divided into segments to boost electric 
potential and maximize power output as illustrated in Fig. A.5| (a)-(b). The 



actually geometry and structure of the TE tube can be engineered to have 
various cross-sectional shapes, longitudinal connections, etc, which will not 
be considered here. The interested reader is referred to [TBI HH1 EHI Then one 
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Figure 3: The maximum dimensionless factor / max (/3) and the corresponding maximizer 

^max(/^) • 



can find that each TE tube has a power output given by (Appendix): 

P f 5T 2 L th l th l-exp(-A) 
° ut " 4t tb ' A 

where i t b is the (wetted) perimeter of the TE tube, A = 7tb ' \ ' t^A^ > 
C p = 4.18 x 10 3 J /kg ■ K is the specific heat of water, Pf = as 2 is the power 
factor of the TE material, n q = k + To<rs 2 /2, and k, a,s is the thermal 
conductivity, electric conductivity, Seebeck coefficient of the TE material, 
respectively. Therefore, in account of the power loss Pi oss (cf., (|2j)), the 
power gain from each TE tube is given by 

12 D„ VT 2 7 



Qi h _ P f 6T% h L tb 

where 



Pgain = Pout " ^TtotW^ = A+ ~ /(A, P), (4) 



,(A,fl = i^3<=*>-£, (5 ) 
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n/ — rv., _|_ 7pp^p P A 2 b Ar tb a _ 16 ^ K g . L tVtb j },„„„„ \ _ 1 . 7tb . 

7tot - 7tb + 7tb L tbA 2 p iv pp ' P - ^ot p2c2pPfST 2 and nence A - s — 

From the above equation it is clear that a large power factor Pf is crucial 
to improve P ga m- The power factor can be significantly improved by com- 
posites. To see this, we consider laminates of two materials as illustrated 
in Fig. [4] (a). From the definition Pf = as 2 , we see that the power factor 
will be improved if electric conductivity increases while Seebeck coefficient 
s is retained. We therefore consider laminates of good conductors, e.g., Cu 
(a = 5.85 el /Vim, k = A01W/mK, s = 1.90fxV/K) or Al (a = 3.66 el /Vim, 
K = 237 W/mK, s = —1.66fiV/K) [9], with TE materials with large Seebeck 
coefficients, e.g., Bi-Te semiconductor systems of p-type doped (Bi .258^.75)2 
Te 3 (a = 0.33 eh/Vim, k = 0.559 W/mK, s = 2±5(jV/K) or n-type doped 
Bi 2 ( Teo9 4 Sboo6)3 (" = 0.36 e5 /ttm, k = 0.788W/mK, s = -209/jV/K) 

Since the electric conductivity and thermal conductivity between the met- 
als and semiconductors differ by more than two orders of magnitude, for the 
applied boundary conditions illustrated in Fig. [5] (a) we anticipate that elec- 
tric potential and temperature drop mainly occurs across the semiconductor 
layer, and hence the Seebeck coefficient remains roughly to be that of the 
semiconductor, i.e., the effective Seebeck coefficient s e = s sia . For uncou- 
pled electric or thermal conductivity problems, it is classic that the effective 
conductivities of laminates are given by 

" =( + ) , K E = (— + ) \ 

"mt ""sm ^mt ^sm 

and hence the effective power factor is given by 

P] = a%sr = (— + —r 1 (s sm y, (6) 

"mt "sm 

where 8 m t and 8 sm are the volume fractions of metal and semiconductor in 
the laminate, respectively. The above formula predicts that the power factor 
would increase significantly if 8 m t increases since a e increases monotonically 
from that of the semiconductor to that of the metal as 8 mt increases from 
to 1. 

Equation ^ is accurate when 9 mt is not close to one. On the other hand, 
it is desirable to have 9 mt close to one to improve power factor as much as 
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Figure 4: Power factors of simple laminates of two materials: (a) a sketch of the laminate, 
and (b) the power factor as a function of volume fraction of the metal. The solid curves 
are plotted according to the self-consistent continuum theory, i.e., the dashed lines 
are plotted according to the rough estimate There is no discernible difference between 
two predictions for lower volume fraction of metal. 



possible. In account of electric and thermal couplings, a continuum theory 
for TE materials has been systematically developed in Liu (2012), where the 
effective TE properties of the laminate are predicted as 



' T a e T$a e s e 
T$a e s e T$[K e + T a e (s e ) 2 }_ 

Tq(T mt Tq 
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Figure Q (b) shows the power factors of Cu and p-type doped ( Bi 2 5Sb .75)2Te 3 , 
Cu and n-type doped Bi 2 ( Teo.94Sbo.o6)3, Al and p-type doped (Bi .25Sbo.75)2Te3, 
Al and n-type doped Bi 2 ( Teo.94Sbo.o6)3 as a function of volume fraction 6 mt . 
We observe that there exists an optimal volume fraction, typically very close 
to one, such that the power factor is maximized with the maximum power fac- 
tor about 100 times larger than the constituent semiconductor. In particular, 
for Cu and p-type doped ( Bi .25Sbo.75)2Te 3 , the sandwich laminate has max- 
imum power factor Pf = 0.206 W/mK 2 (a e = 8.30e6 /Sim, n e = 158 W/mK, 
s e = 145 fiV/K) when #cu = 0.9975. We also remark that the sandwich 
structure illustrated in Fig. [41 a) is convenient for mass-production and for 
electric connections since the metal can serve as electrodes as well. 

Further, from Q it is clear that small f3 is favored for large /(A, (3). We 
shall assume (3 < 1 by tuning geometric factor L t b, hb, A;t>- Upon maximizing 
f(\,(3) over A > for given (3, we obtain the maximum f ma , x (f3) and maxi- 
mizer A max (/3), which are shown in Fig. [3] (c). From Fig. [3] (c) we observe that 
/max depends on (3 weakly, varies from 1 to 0.1 when (3 varies from 10~ 5 to 
1. Moreover, by Q the power gain from each TE tube P ga in scales as l/t t b- 
Therefore, a thin tube wall is preferred and the smallest ttb is limited by the 
manufacture and structure requirement which may be reasonably chosen as 
t t b — 0.1 ~ 1 mm. 

To fix the design parameters, we select as the TE tube wall material 
the sandwich laminate of Cu and p-type doped (Bi .25Sb .75)2Te 3 with the 
maximum power factor Pf. For maximum power factor, the thickness of 
semiconductor layer between metals shall be chosen as t sm = t t b * 0.25% = 
0.25 ~ 2.5 /im, which does not appear to be a challenge for modern technol- 
ogy of nanofabrication. Assume that 5T = 10K which suffices to account for 
the heating up of cold water when being pumped up from deep sea. Then 
the dimensionless parameters (3, A are given by (SI units): 

^tb^tb 7tot h 

Assume both the long underwater pipe and short TE tubes are circular and 
select i?pp = 1 m, L pp = 10 3 m, iVpp = 1, i? t b = 10 -2 m, L t b = 1 m, 
A t b < 10 5 , and hence 7 tot ~ 7tb = 87r, (3 = 2.21 x 10~ 5 /t tb . Upon choosing h 
such that A = A max (/3), we find the power gain P ga i n = 0.324 • / max /t t b Watt 
from each TE tube and the number of TE tubes needed for a megawatt 
power gain which are listed in Table [T] We remark that rectangular fin-like 
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TE tube may be more favorable in practice for ease of fabrication and mass- 
production which, however, has little effect on our prediction of power gain 
from each TE tube. 

Designs of a megawatt TE generator (ST = 10K). 

Physical means and dimensions of the symbols: ttbj Rtb, L t b — thickness, inner radius 
and length of TE tube in mm, cm, m, respectively; / max the maximum dimcnsionless 

factor in Q; h in meter — the pressure needed to drive the flow in pgh (cf., Fig. [I]); 
-fgain — the power gain from a tube in KW; N t b — the number of tubes needed for one 

megawatt power gain. 

Table 1: Proposed designs of a megawatt TE generator (ST = 10K) 2 . 
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3. Conclusion 

The estimated cost per year of one megawatt electricity in 2016 is about 
0.83 million USD from conventional coal plants to 1.84 million USD from 
photovoltic power plants |15j . Also, the investment of the world largest pho- 
tovoltic farm (Sarnia PV power plant in Canada (22]) is 300 million USD 
for an annual capacity 1.2 xlO 5 MWh, or equivalently, an effective power of 
14 megawatt (designed peak power 60 megawatt). Even if Table [l] under- 
estimates the number N pp of TE tubes by 10 fold, it still appears feasible 
to manufacture a megawatt TE generator using waters of 10K temperature 
difference as "fuel" which lasts for 20 years and costs less than 20 x 1.84 mil- 
lion USD, and is hence economically competitive compared with photovoltic 
farms. For geothermal sources, we anticipate the temperature ST is at the 
order of 50 K or higher, implying an increase of power gain by 25 fold since 
-Pgain scales as ST 2 . 

In summary, we have presented a simple design of heat exchanger or elec- 
tricity generator which use heat resources of small temperature differences. 
The theoretical analysis shows that it is feasible and economically competitive 
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to build large-scale ocean-based or geothermal power plants by thermoelec- 
tric effects. A for-the-moment science-fiction vision is that power plant units 
like barges are manufactured on shore and towed to a selected tropical ocean 
site where electricity is generated and transmitted to shore. Whenever more 
power is needed, new units can be manufactured, towed and assembled to 
the site. The renewability and scalability of the TE power plants may yield 
a way to resolve the energy crisis facing human society. 



Appendix A. Flow and thermoelectric problems in a TE tube 

Appendix A.l. Flow problem 

For a low Reynold number (Re < 2300) steady flow in a circular pipe of 
length L and radius R and maintained by a pressure drop AP, the flow is 
laminar and referred to as the Hagen-Poiseuille flow. Let v = v(r) be the 
velocity of the flow, and by the Navier-Stokes equation we have 

v(r) = if r = R, 

where \i is the viscosity of water at 293K and p' = —AP/L is the pressure 
gradient. The solution to the above problem is given by 

, , AP(R 2 - r 2 ) 
V{r) = A/iL ' 

and henceforth volumetric flow rate is given by 

( K M „ , 71 APR 4 
Q= / v(r)2ixrdr = — -— . (A.l) 
Jo 

The power dissipated by the viscosity is given by 

R ,„ ,2 , , 7lAP 2 R 4 

[J>\vv\ znrdrdx = 



^0 



8/iL 



)2 



= APg = 8j ^~ =: VlLQ 2 /A 2 , (A.2) 
where A is the cross-sectional area and 7 = 811 for circular pipe. 
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Appendix A. 2. Thermoelectric problem 

We consider the thermoelectric problem as water flows through the TE 
tube immersed in an ambient medium of 10 K temperature difference. The 
TE tube will be divided into multiple segaments to boost electric poten- 
tial difference and power output. Assume that each segment of tube is short 
enough such that the temperature and electric potential on interior and exte- 



rior walls may be regarded as constants on each segment, see Fig. A. 5 a)-(b). 
Further, it is desirable that each segment of the TE tube works at the opti- 
mal conditions to have maximum power generation per unit volume, i.e., the 
segment length U and total electric current J shall be controlled such that 
the electric current density across the wall of each segment is given by (cf., 
§ 3.1.1 of Liu, 2012) 

= ggkpjM , (A.3) 

where the temperature of the exterior walls of all TE tubes is assumed to 
be constant and equal to T ex , and X{ is the midpoint of the iih segment. 
Since the electric current in the closed circuit is constant, the length of each 
segment U shall be proportional to 1/(T CX — T(xi)). At this condition, the 
efficiency of conversion is not optimal and given by 

V = ^carnot, Vs = 4 + 2ZT ' 

which is only 0.5% less than the optimal efficiency rj° pt = (\/ZTq + 1 — 
l)/(\/ZT + 1 + 1) for ZT = 1, and even smaller for smaller ZT . Also, the 
potential difference between exterior and interior walls at the ith segment is 
given by 

A*ex(^i) -£*in(zi) = -|(T ex -T(Xi)) (i = l,--- ,N). (A.5) 

The reader is invited to check for the detailed calculations of above results 
in Liu (2012 [23]). 

Moreover, the heat flux flowing out of the tube across the tube wall is 
given by 

q(x) = - Kq Tex ~ T{X \ K q = K + T,as 2 /2. (A.6) 

ttb 
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(cl ) Electric connection ofthe TE tube of the same type 
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Figure A. 5: (a) Electric connection of the TE tube in the longitudinal direction if the 
type of TE material remains the same; (b) Electric connection of the TE tube in the 
longitudinal direction for alternating type of if TE materials. 



For simplicity we assume that the length of the tube is much larger than the 
cross-sectional diameter of the tube, and hence temperature may be approx- 
imately regarded as uniform on each cross-section of the tube. Neglecting 
the heat generated by viscosity, by the energy balance we have 



pCpQtb-^T^x) + l tb q = 0, 



(A.7) 



where Z t b is the length of the perimeter of a cross-section of the tube. In 
account of the inlet temperature T(x = 0) = T° which may be assumed to 



be the temperature of the source, by (A. 6) and (A.7) we have 



T(x) = T ex + (T° - T ex ) exp(-x/L* 



(A.J 



where 



pCpQtbttb 
Kqhb 



(A.9) 



13 



Therefore, the outlet temperature is given by 




(A. 10) 



L* pCpQtbttb 



The power generated by each tube is a summation of power generated by all 
segments and well approximated by (cf., § 3.1.1 in Liu, 2012) 
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